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The most massive galaxies in the present-day Universe are found to lie in the 
centres of rich clusters. They have old, coeval stellar populations suggesting 
that the bulk of their stars must have formed at early epochs in spectacular 
starbursts 1 — luminous phenomena at submillimetre wavelengths 2 . The most 
popular model of galaxy formation predicts that these galaxies form in proto- 
clusters at high-density peaks in the early Universe 3 . Such peaks are signposted 
by massive high-redshift radio galaxies 4 . Here we report deep submillimetre 
mapping of seven high-redshift radio galaxies and their environments. These 
data confirm not only the presence of spatially extended massive star-formation 
activity in the radio galaxies themselves, but also in companion objects previ- 
ously undetected at any wavelength. The prevalence, orientation, and inferred 
masses of these submillimetre companion galaxies suggest that we are witness- 
ing the synchronous formation of the most luminous elliptical galaxies found 
today at the centres of rich galaxy clusters. 
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Whilst existing submillimetre studies of high-redshift radio galaxies (hereafter HzRGs) 
have shown that their star-formation rates are large enough to build a massive galaxy in 
< 1 Gyr 2 ' 5 ' 6 ' they have provided no information on the spatial extent of this emission 
or on the star-formation activity in their environments. We have therefore mapped a 
sample of seven objects with redshifts ranging from 2.2 to 4.3 at a wavelength of 850 jum 
with the Submillimetre Common-User Bolometer Array (SCUBA) 8 on the James Clerk 
Maxwell Telescope (JCMT). The targets were chosen from those sources found to be 
submillimetre bright in the previous SCUBA surveys of HzRGs 6 ' 7 . Our new maps illustrate 
the distribution of dust-reradiated emission in and around the HzRGs on scales from 5" 
to 160", or 30kpc to 1 Mpc. We illustrate the seven submillimetre maps from this survey 
in Figure 1; the orientation of the radio jets of each HzRG is represented by tick marks on 
these maps. 

One of the most striking aspects of the submillimetre maps is that the dust emission 
from the central radio galaxy is resolved in at least five of the seven sources - even with 
the coarse beam of the JCMT. In Figure 2 and Table 1 we present simple two-dimensional 
Gaussian fits to the data which, while not giving a true reflection of the physical situation, 
at least provide a quantitative measure of the spatial extent of the dust emission. This 
emission is sometimes in the form of several partially-resolved or merged clumps (typified 
by 8C 1909+722), sometimes in an apparently smoother distribution (e.g. 4C 60.07), and 
is more extended than the radio emission in most cases. The extent of the dust emission 
ranges from 50 to 250 kpc, a physically interesting size because (1) the corresponding 
half-light radii (30-150 kpc) are equivalent to those of brightest cluster galaxies in the 
local Universe 9 , and (2) gas-dynamical simulations of major galaxy mergers predict that 
star formation should peak when the galaxies are separated on approximately this scale. 10 
Higher resolution, high signal-to-noise millimetre/submillimetre imaging observations are 
required to investigate these possibilities. 

A second striking feature of our SCUBA maps is that several of the fields also contain 
serendipitous detections of new, luminous submillimetre galaxies. In Table 1 we list the 
properties of those companion galaxies with peak signal-to-noise ratio > 4. Note that 
several of these companions are, like the HzRGs, resolved at 850 ^m. The source density 
of companions is also higher than found in 'blank-field' surveys; current submillimetre 
number counts 11 indicate that we should find, on average, about one random submillimetre 
source per SCUBA field at a flux level of Ssso > 5-6 mJy compared to our detection of 
approximately twice this number (Table 1). Furthermore, the 4C 41.17 and 8C 1909+722 
fields each contain a companion source with Ss^o > 20 mJy. By comparison, the biggest 
blank field survey conducted to date 12 , which has an area ~ 7 times that of our own, 
contains no robust sources with Ss5o > 15 mJy. 

Arguably the most surprising feature of the maps is the observed aligments between 
the HzRG radio axis and (1) the submillimetre emission (2) the brightest submillimetre 
'companions'. The significance of these alignments is presented and quantified in Figure 3. 
These results can be interpreted as follows. First, while one might be tempted to conclude 
that the first alignment effect is indicative of jet-induced star formation similar to the 
radio-UV alignment effect previously reported in HzRGs 13 ' 14 , the fact that in all but two 
cases, the radio source is also aligned with submillimetre source on scales well beyond the 
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size of the radio emission indicates that this is unlikely to be the complete explanation. 
Moreover, this model clearly cannot explain why the sub millimetre emission from the 
HzRGs appears to be aligned with positions of the brightest companion sources in their 
vicinity (in at least 4 out of 7 cases). 

There is one model that can explain both alignment effects. We propose that the 
brightest submillimetre companions trace the large-scale structure around the HzRG and 
that these directions thus contain the densest cross-sections of gas. Next, by selecting some 
of the very brightest known radio sources at this epoch we will then have (inadvertently, 
but probably inevitably) selected sources which happen to have produced jets aligned 
with the densest regions of gas, thereby producing very effective working surfaces and the 
brightest hot spots. Such a selection effect has been suggested before to explain apparent 
large-scale optical-radio alignment effects at lower redshift. 15 ' 4 

Of course, an important corollary to this interpretation is that the brightest com- 
panion sources seen in the submillimetre maps must lie at the same redshift as the radio 
galaxies, occupying the same large-scale structure. However, the fact that extension of this 
analysis to the second brightest apparent companion object does not reveal a significant 
alignment effect shows that not all of the objects seen in these images need lie at the same 
redshift as, and be physically associated with the radio galaxies. This is consistent with 
our earlier estimate of the rate of contamination in these fields by unrelated submillimetre 
sources. Thus the analysis of these possible alignment effects indicates that the submillime- 
tre emission from the radio galaxies themselves, and that of the brightest companions, is 
tracing the large scale structure around radio galaxies at z ~ 3. We note that for 53W002, 
a HzRG not included in our sample, it has been shown with optical spectroscopy that the 
brightest submillimetre companion is indeed at the same redshift as the radio galaxy. 16 A 
similar conclusion can be inferred for PKS 1138 — 262 where the brightest submillimetre 
companion is coincident with one of 5 actuve galactic nuclei (AGN) forming a radio-aligned 
filament in the plane of the sky, and which has the same redshift as the HzRG based on 
narrow-band imaging of redshifted Lyct emission. 17 

Since we have shown that these fields contain overdensities of submillimetre sources, 
it is of interest to ask what is the typical mass of these galaxies. We can do this by 
assuming that they lie at the same redshift as the HzRGs. The dust masses of the HzRGs, 
calculated from their 850- jum fluxes 18 are given in Table 1. Using standard assumptions 
we can convert these dust masses into total gas masses, and hence estimate the masses of 
their associated dark matter haloes (see table caption for details). We estimate that they 
reside in dark matter halos with masses in excess of 10 12 M . Since the submillimetre flux 
densities of the companions are similar to those of the HzRGs we can infer that they also 
reside in dark matter halos with masses in excess of 10 12 M Q . We therefore conclude that 
these regions can contain more than one very massive and gas-rich galaxy of > 1O 12 M0. 
Can we reconcile this result with the hierarchical models of galaxy and structure formation? 

By design we have selected highly-biased regions of the high-redshift universe by imag- 
ing around some of the most luminous radio galaxies at these epochs, which are expected 
to host some of the most massive black holes. We can attempt to quantify the impact 
of this bias by exploring the predictions of numerical simulations of the collapse of Cold 
Dark Matter haloes within a A-dominated universe. The details and results of these sim- 
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ulations are given in Figure 4. These simulations do indeed appear to be in accord with 
the observations, i.e. both data and theory suggest that the submillimetre companions 
revealed by our SCUBA imaging are associated with dark-matter halos more massive than 
10 12 M©. This conclusion is, however, subject to an assumption about the duty cycles of 
the submillimetre luminous phase; if true it would imply that the vast majority of such 
massive halos are actively forming stars at the epochs sampled by these images. This is 
not unreasonable, especially since we have deliberately targetted fields in which at least 
one massive object (i.e. the HzRG) is actively engaged in intense star formation. 

Pursuing this comparison one step further, we can explore the predicted and observed 
distribution of companion objects in this high- mass range. In fact, while the average 
number of such companions in the simulated SCUBA images centred on the HzRG is 1, 
this average arises from a skewed distribution. As shown in Figure 4, of the 40 high-mass 
haloes investigated, 50% have no such massive companions, with the average of 1 resulting 
from the fact that, of the remaining half, 40% have 2 companions or more. We have 
over-plotted in Figure 4 the corresponding histogram of companion incidence for the 7 
HzRGs imaged in this study (after statistically correcting for field contamination). This 
comparison at least illustrates that the distribution of companion incidence in our data is 
consistent with the prediction of the simulations, i.e. typically we see either no companion, 
or > 2 companions. 

The average number of companions, the distribution of companion incidence, and the 
inferred baryonic gas masses of the companions are all most consistent with the interpre- 
tation that the SCUBA sources uncovered in this study are the progenitors of massive 
present-day cluster ellipticals. In this case, since elliptical galaxies are known to contain 
massive black holes with mass proportional to that of the spheroid, 19 it is reasonable to as- 
sume that the black hole and stellar mass grow coevally from the same gas reservoir. This 
in turn suggests that the companion objects should contain buried AGN. 20 ' 21 First results 
show a high rate of correspondence between the submillimetre companions and luminous 
X-ray sources, 16 ' 22 suggesting that this is indeed the case. We are presently pursuing high- 
resolution follow-up observations of these fields at millimetre and optical/near-infrared 
wavelengths. The former will yield information on source structure while the latter will 
reveal counterparts for spectroscopic redshift determination on 10-m class telescopes. 
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Tables. 

Table 1. Target fields. 



Target 


n a- t a 
Coordinates 


(J2000) 


/J 


ab 

peak 


ab 

°total 


M d 


d 

FWHM/P. A. 


Radio size/P.A 


name 


R. A. 


Dec. 




(mjy) 


(mjy) 


(MQ) 


/'A 

( / ) 


( / ) 


8C 1909+722 


19 08 23.30 


+ 72 20 10.4 


3.54 


20.5il.2 


34.9lh3.0 


1.2X10 9 


26X9Q18 


14Q15 


1. 


19 08 27.47 


+ 72 19 28.0 




17.6ztl.2 


23.0^12.5 




12 X 8Q169 




2. 


19 08 29.31 


+ 72 20 49.6 




6.5ztl.2 


8.71L2.4 








3. 


19 08 16.12 


+ 72 20 24.0 




6.5ltl.2 


4.3lb2.1 








8C 1435+635 


14 36 37.33 


+63 19 13.1 


4.261 


5.9ztl.O 


6.oi2.1 


3.3X 10 


7Xo@159 


4Q156 


1. 


14 36 32.46 


+63 20 02.5 




4.3lLl.O 


3.9ztl.8 


i.oXio 9 






B3 J2330+3927 


23 30 24.91 


+ 39 27 11.2 


3.086 


15.6il.l 


22.2lh2.7 


16X7Q130 


2Q149 


1. 


23 30 19.14 


+39 27 03.0 




6.9il.l 


8.2ltl.9 








2. 


23 30 20.52 


+39 26 57.7 




s.oii.i 


3.5ztl.4 








PKS 1138-262 


11 40 48.25 


-26 29 10.1 


2.156 


5.9il.l 


6.71L2.4 


4.6 X 10 


30Xo@72 


111986 


1. 


11 40 53.38 


-26 29 11.9 


~2.16 


9.9ztl.l 


7.8lb2.2 




17X0Q60 




2 




or on Ka a 




5 Q — 1 — 1 1 


3 1 ~l~ 1 6 








3. 


11 40 45.61 


-26 29 06.6 




4.6lLl.l 


2.2lLl.4 








4C 60.07 


05 12 54.80 


+60 30 51.7 


3.788 


21.6il.3 


23.8lh3.5 


1.3X10 9 


11 X 4Q110 


16O107 


1. 


05 12 46.52 


+60 30 35.8 




5.9ztl.3 


6.3zt2.1 








4C 41.17 


06 50 52.15 


+41 30 30.8 


3.792 


12.3ztl.2 


12.oi2.3 


7.2X10 8 


9X6897 


20Q48 


1. 


06 50 51.52 


+41 30 01.5 




12.2ztl.2 


21.2lt2.9 




17X 15Q54 




2. 


06 50 49.25 


+41 30 01.5 




7.1ZC1.2 


6.2lbl.9 




8 X 0Q50 




WN J0305+3525 


03 05 47.42 


+35 25 13.4 


3±1 


12.2311.0 


17.6lL2.4 


6.9X 10 


19 X 9@9 


2Q64 



a Coordinates of the HzRG companion sources are measured from the submillimetre maps. They are 
accurate to 4—5 arcseconds. 

b (Speak is tne nux density per beam calculated at the peak of the dust emission. Stotal is the flux density 
of total dust emission calculated in an aperture. A calibration uncertainty of ~ 10% is not included 
in these values. For some sources, the higher total flux densities reflect the extended nature of the dust 
emission, consistent with the finding that interferometric millimetre flux densities are often smaller than 
those calculated from single dish measurements. ^8>23 

c Dust mass calculated from the peak fluxes, and assuming a dust emissivity index, /3=2.0; a dust tem- 

o 1 

perature, T c j=40 K typical of HzRGs; a mass absorption coefficient, k(j(850 /Um)=0.076 m kg ; and an 
=0.3, Qa =0.7, h =0.7 cosmology. Our selection criteria result in a small spread in the calculated 
dust masses; an average value is ~ 8 X 10^ Mq. We can convert this dust mass into a conservative 
gas mass (molecular and atomic) assuming a standard Galactic gas-to-dust mass ratio of 200 2 ^ giving 
1.6 X 10 11 M . We can confirm this estimate using the measured CO luminosities of three of the HzRGs, 
which can then be turned into molecular gas, and then total gas, mass using the standard Galactic con- 
version factor between CO luminosity and H2 mass and a ratio of atomic-to-molecular gas of ~ 2. In 
this manner, gas masses of 1.3 — 2.2 X 10 11 M have been estimated for 4C 60.07, 8C 1909+722 and 
B3 J2330+3927. Hence the gas masses from the two methods are in reasonable agreement, and a 

representative estimate for the total gas mass in a HzRG with a submillimetre flux density of 10-15 mjy 
is thus ~ 2 X 10 11 Mq. While this number is arguably uncertain by a factor of a few, it is very hard to 
argue that it can be inflated by an order of magnitude. 5 Assuming an order of magnitude more mass is in 
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the form of dark matter we estimate the halo masses to be > 10 Mq. 

d From a 2-dimensional Gaussian fit to the 850 /im emission, deconvolved from the beam. All P.A.s are 
measured east of north. For the companion galaxies we give sizes for the brightest (highest S/N) sources 
only. A quoted size of zero arcseconds means that the source was not resolved in that direction. 
e Quoted radio sizes are hotspot separations measured at 5 GHz. 

f WNJ 0305+35 lacks a robust spectroscopic redshift, but is thought to lie at Z = 3 i l. 30 
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Figure Captions. 

Figure 1. Continuum emission from dust in and around seven high-redshift radio galaxies. The data for 

26 8 
4C 41.17 have been published previously. Data were obtained using the SCUBA submillimetre camera 

on the James Clerk Maxwell Telescope (JCMT) in its dual-wavelength mapping mode during 1998-2001 

in the top 20 percentile of Mauna Kea weather conditions. The beam, after smoothing to a FWHM of 

14.4" at 850 /Um, is shown bottom right in the form of a map of the blazar, 3C 345, with contours at —40, 

—30, —10, +10, +30, +50, +70, +90% of the central peak. The JCMT secondary mirror was chopped 

and nodded by 30" in Right Ascension (R.A.), resulting in the — 0.5/+1.0/ — 0.5 beam profile. After all 

overheads, 30-40 ks were spent integrating on each of the fields, a total of ~100 hr before overheads, split 

down into 1280- s integrations separated by checks on focus, pointing accuracy and atmospheric opacity. 

All images have been cleaned with the smoothed beam. Submillimetre contours are shown at —3, 3, 4, 

5, 6, 8, 10 X(T, where <7 includes the contribution due to sources below the detection threshold and to the 

chopping and nodding procedure and ranges from 1 to 2 mjy beam - . The galaxy with the most compact 

appearance, 4C 41.17, actually subtends ~7.5" FWHM when deconvolved from the beam. The tick marks 

around the central radio galaxies show the direction of the kiloparsec scale radio jets. 

Figure 2. Submillimetre extent of the radio galaxies. Radial profiles of the dust emission compared 
with that of the beam (solid line). The x-axis shows the radius from the centre of the radio galaxy dust 
emission (arcseconds) and the y-axis shows the normalised flux density. An accurate beam profile or point 
spread function (PSF) was measured with a ~ 1 hr mapping observation of the blazar 3C 345 which 
had an 850/Xm flux density of > 2 Jy at the time. The image (shown in Figure 1) has a full width 
half maximum (FWHM) of 14.4". All subsequent size measurements are deconvolved from the PSF (see 
Table 1). The reality of the extended nature of the dust emission is established from observations of 
the quasar, BR 1202 — 0725, using the same technique and integration times used for the radio galaxies 

to mimic any systematic errors. BR 1202 — 0725 serves as a particularly conservative comparison object 

27 // 
— its dust continuum was resolved by the IRAM interferometer into two components separated by 4 

(P.A. 120il0°). As expected, it appears relatively compact (6.5" X 5.5" at P.A. 6°) in our submillimetre 

maps, and its radial profile is much more similar to that of the beam than to those of the HzRGs. We 

conclude that our observing procedure does not produce spurious morphologies, that the extended nature 

of the submillimetre emission from the radio galaxies is real and occurs on scales that are typically 3>4". 

Figure 3. Observed alignment effects. The histograms show the offset in degrees between the kilo- 
parsec scale radio structure of the HzRG and a) the position angle of the submillimetre source, b) the 
position angle defined by the projected vector joining the radio source to the brightest, spatially distinct 
submillimetre source, and c) the position angle defined by the projected vector joining the radio source 
to the second brightest, spatially distinct submillimetre source. In four cases these submillimetre sources 
are the companions listed in Table 1. In the remaining three cases one or both of the companions has a 
signal-to-noise ratio of 2 to 4; while these sources are not listed in Table 1 we use them in the analysis for 
completeness. We have applied the Kolmogorov Smirnov test to assess whether the apparent alignment 
effect seen in the first two of these plots is indeed statistically significant. Despite the small number statis- 
tics, we find that both of these apparent alignment effects are significant at the 2-<7 level (as compared 
to a random distribution; p = 0.04 and p = 0.05 respectively) whereas the corresponding result for the 
third distribution is not significant (p = 0.3). 

Figure 4. Comparison with hierarchical models. The predicted distribution of number of companion 
sources per SCUBA field (after subtracting a background of 1 source per map for field contamination) from 
our A-CDM simulations shown against the observed distribution (shaded region). For the simulations we 
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assume that the radio emission from the HzRG is powered by black holes with mass — 10 9 M© 24 , and 
therefore these lie at the heart of dark matter haloes with masses of ~ 10 13 Mq^. Allowing for intrinsic 
scatter and uncertainties in this estimate, we therefore searched for virialized dark matter haloes with 
masses in the range 12.5 < log^Q^-Mhaloj/M© < 13.5 at Z — 3 within an N-body simulation of a 
region ofaA-CDM universe, covering a comoving volume of (lOO/l 1 Mpc) 3 . Parameters of the simulation 
were matched to the "concordance model" with m = 0.3, Qa = 0.7, h = 0.7; further details of the 
simulation are given elsewhere (ACDMiooa) 2 ^- Within this simulation we identified 40 haloes which 
could be the counterparts of the HzRGs, and then investigated the mass distribution of other virialized 
dark-matter haloes within a sphere of projected diameter 3-arcmin (approximately the SCUBA field of 
view) around each of these. Within these regions we found typically 10 haloes with M > 10 11 Mq, but 
on average only 1 halo with M > 10 12 M . 
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